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When we calculate gravitational waveforms from extreme-mass-ratio inspirals (EMRIs) by metric 
perturbation, it is a common strategy to use the adiabatic approximation. Under that approxima- 
tion, we first calculate the linear metric perturbation induced by geodesies orbiting a black hole, 
then we calculate the adiabatic evolution of the parameters of geodesies due to the radiation reaction 
effect through the calculation of the self-force. This procedure is considered to be reasonable, how- 
ever, there is no direct proof that it can actually produce the correct waveform we would observe. In 
this paper, we study the formal expression of the second order metric perturbation and show that it 
be expressed as the linear metric perturbation modulated by the adiabatic evolution of the geodesic. 
This evidence supports the assumption that the adiabatic approximation can produce the correct 
waveform, and that the adiabatic expansion we propose in Ref.[l|] is an appropriate perturbation 
expansion for studying the radiation reaction effect on the gravitational waveform. 
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I. INTRODUCTION 



We study the gravitational waveforms from particles moving around Kerr black holes by using a metric perturbation 
method. There is an established method for calculating these waveforms from the linear metric perturbation of 
1 q' 1 Kerr black holes. By the consistency of the Einstein equation, the source stress-energy tensor of the linear metric 
perturbation must satisfy the conservation law with respect to the background. As a result, the source of the linear 
metric perturbation moves along a geodesic of the background Kerr metric. Due to the integrability of the geodesic 
(3J[)! equation for the Kerr metric, its bound solutions are stable and have periodic features. To understand some feature 
of gravitational waves from such a stable system, we develop a technique of formal calculation 0. The advantage of 
| this technique is that, one can easily grasp some key features of the waveform without a complicated calculation, and 
it helps us to construct the strategy for an explicit calculation. Using this technique, we find that the waveform from 
' the linear metric perturbation has a periodic feature[2[ as we review in SecfllJ Based on our understanding of this 
feature, a present numerical code is trying to identify which of the gravitational wave modes are strong enough to be 
observed by gravitational wave detectors [3] . 

Because gravitational waves carry away energy and angular momentum, the system must have a dissipative evo- 
lution. This should change the periodic features of the gravitational waveform that we know from the linear metric 
perturbation. New features must be seen in the second order metric perturbation because the source term must 
include this dissipative effect. The goal of this paper is to find these new features. That is, the goal is to find how the 
wave amplitude and phase are modulated by gravitational radiation reaction. However, a serious calculation of the 
second order metric perturbation still has a lot of technical difficulties. We therefore extend the technique of formal 
calculation we developed to study the linear metric perturbation. We expect that the result of this formal calculation 
will be useful in future explicit calculations of the second order metric perturbation. 

This problem has attracted the attention of the gravity community[J| because its relevance to gravitational wave 
detectors, especially LISA. Among the primary targets for LISA are extreme-mass ratio inspirals (EMRIs), the 
inspiralling binary systems of supermassive black holes (with mass ~ 10 5 — 10 9 Mq) and stellar mass compact objects 
(with mass ~ 1 — IOOMq). Because of the extreme mass ratios, the metric perturbation of the black hole is effective 
for studying the dynamics of the system; we use the Kerr geometry of the supermassive black hole as a background, 
and approximate the compact object as a point particle. The linear metric perturbation may predict the waveform 
at an instant. However, LISA will detect gravitational waves for several years. This time scale is comparable to the 
radiation reaction time scale of EMRIs and is essential to consider the dissipation effect of the gravitational radiation 
reaction because the effect will accumulate during the observation time. 

Most of the relevant investigations made so far discuss the calculation of the self- force [f| which needs regularization 
Q. The underlying idea for such calculations is that: (1) At each instant, the orbit is approximated by a geodesic 
since, due to the extreme mass ratio, gravitational radiation reaction is a small effect. (2) As the effect of radiation 
reaction accumulates over time, the orbit changes from one geodesic to another. (3) These changes can be deduced 
from the self-force. This idea is usually referred as the 'adiabatic' approximation. This calculation strategy has some 
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theoretical problems. Since the self-force is gauge dependent, we may not be able to make a unique prediction about 
the orbital evolution. In Ref. [l[, we showed that, for times over which the standard metric perturbation expansion 
is valid, the self-force can be arbitrarily adjusted via gauge changes. As an extreme example, we showed that there is 
a gauge transformation which completely eliminates the self-force. 

Do we still need to calculate the self- force? In Ref. we argue that the answer may be yes. Under a certain gauge 
condition, the self-force may include the correct radiation reaction effect With this gauge condition, the calculation 
by the 'adiabatic' approximation will give us the correct prediction of the gravitational waveform. The present paper 
discusses the self-force under that 'physically reasonable' class of gauge conditions proposed in Ref. [l], 0] in the 
context of the second order metric perturbation. We will focus on gauge invariant quantities. The basic idea is that, 
because the gravitational waveform is observable, the features we can see from the waveform have invariant meanings. 
We will show that the self-force correctly describes the radiation reaction effect through the second order metric 
perturbation. 

The field equations for the linear and second order metric perturbations are derived as follows. We expand the 
metric in a small parameter e 

9^ = 9$ + eg$ + e 2 9$ + 0(e 3 ) , (1.1) 

where gffi is the background metric. We insert the e-expansion of the metric into the Einstein tensor and expand it 
in powers of e. We formally obtain 

qv* = eGW^^Wj+e^GW^fe^l+G^tgW^WlJ+OCe 3 ), (1.2) 

where G^^h] is linear in /i M „ and G^^hW, h' 2 ^] is bi-linear in hP~^ v and hffi. By similarity expanding the source 
term 

Tiiv = eT [l]F + e 2 T [V + ( e 3) ] ( 13 ) 

we can formally write the perturbed Einstein equations to second order e 

G [i]^[ g (i)] = SttTW^, (1.4) 

G [1] Hg (2) ] = -G [2] Hg (1) ,g (1) ] +SnT^. (1.5) 

In Secim we review the derivation of ljl.4|) [2j. Due to the consistency of the Einstein equation, the source of the 
linear metric perturbation, T^^ v ', must be a geodesic. Because we are interested in a bound geodesic, bound geodesies 
are triperiodic, and we will see that the linear metric perturbation induced by it is as well. Before discussing the 
second order equation (|1.5|) . we make a formal argument about the self- force and the orbital evolution due to this 
effect in Sec. Mil because those effects should be included in the source of the second order equation, T^^. In Sec. 
IIVI we discuss the formal analysis of the second order Einstein equation (jl.5p . The second order metric perturbation is 
induced by the quadratic term of the linear metric perturbation, and by the second order source term. We emphasize 
these two effects leads to substantially different features of the second order metric perturbation. Sec. [V] concludes 
our result. 

Our background geometry will be a Kerr black hole with mass M and spin parameter a. We will a point particle 
with mass [i as the source of the metric perturbation. Throughout this paper, we use Boyer-Lindquist coordinates, 
{t, r, 8, <fi} — {x a }. We adopt the geometrized units which are defined such that G = c = 1. 

II. LINEAR METRIC PERTURBATION 

In this section, we review the periodic feature of the linear metric perturbation. It is well known that, when the 
background is a vacuum solution, the source term of the linear metric perturbation must be conserved with respect 
to the background. Since we use a point mass source, this means that its world line is a geodesic of the background 
geometry. Because we are interested in gravitational waves from binaries in the inspiralling phase, we only consider 
bound geodesies. 

We denote the orbital coordinates of the geodesic by {x a } and the 4- velocity by v a := dx a /dr, where r is the proper 
time. The geodesic equation for a Kerr black hole has three nontrivial constants of motion: the energy E — firj^v 01 , 
the z-component of angular momentum L — iiv^v a and the Carter constant G = {^/2)rj a pv a v 13 . Here rj^ and rfc are 
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temporal and rotational Killing vectors, respectively, and r\ a p is the Killing tensor. Hereafter we collectively denote 
these constants by £ a . The geodesic equations can then be formally written as 

" lrV R(f ',£*), =e(6;£ a ), (2.1) 



[dXj " v '~ " \dX) 

Jl = T r (f;£ a ) + T e (e;£ a ), g = $ r (f;£°) + ®e{9;£ a ) , (2.2) 

where A is the affine parameter of the orbit related with the proper time by A = Jdr/(f 2 + a 2 cos 2 #). For bound 
geodesies, the r-motion and the ^-motion are periodic and one can write them as a discrete Fourier series. We formally 
write 



r = $> fc e ifc *", Xr = T r (A-A r ), (2.3) 

k 

9 = 5>e«™, Xe = T e {X-X e ), (2.4) 
i 

where A r and X e are constants for integration. The expansion coefficients, and Q\, and the effective frequencies, T r 
and Yg, are functions of the constants of motion £ a . Using (|2.3|) and (|2.4p . we can formally integrate ()2.2jl as 

i = Xt + J2 + E ^ > M = T '( A - A *) > ( 2 - 5 ) 

k I 

^ = X4> + E + E $ eUXe - » = T *( A " ' ( 2 - 6 ) 

fc z 

where A* and A^ are constants of integration. The expansion coefficients, t e k: cff k and cf)f , and the effective frequencies, 
T t and T^, are functions of the constants of motion £ a . In summary, one can specify a geodesic by three constants 
of motion £ a and four constants of integration X a . Because one is free to choose the zero of the affine parameter A, 
one of four integral constants A Q is not physically significant, and the gravitational waveform induced by the geodesic 
depends only on three differences of the four integral constants [See (|2.14j) .] . 

Under a certain gauge condition, one can define a tensor Green's function for the linear Einstein equation (|f .4[) . 
Because the background is stationary and axisymmetric, this tensor Green's function has the form 

G aPliV (x, x') = G al3 ^(t -t',<p- 4>'\ r, r', 9, 6') . (2.7) 

For the linear metric perturbation, the source term is described by the stress-energy tensor of a point particle [f| 
written as 

TW^ = /x / drv^ S{x Z^ -, (2.8) 

where \g^ \ is the determinant of the background metric. The corresponding linear metric perturbation is given by 

gW{x) = 8tt M j dX (^\ Gafinv (x, x(X)) v»v v . (2.9) 

Using the formal expression for the geodesic, (|2.3[) - (|2.6|) , we obtain the following formal expression for this linear 
metric perturbation 

&*) = ^ e --(M,-)t+»^I i ; ()[Ara) ( r)fl ) e i(Hf +l«,t ( +™/) ; ( 2.io) 

k,l,m 

r^-i 

U(k,i, m ) = ku; r + lug + mu<j, , u> r = ^ , uj 9 = — ^ , = — ^ , (2.11) 

It I t It 

f = T t (X r — A*) , t e = T t (X e -X t ), i* = T t (A* - A*) , (2.12) 

where the expansion coefficients, <7^"L fe ; m ){ r i 9), depend only on the constants of motion, £ a , and are independent of 
A r , A e , A* and A*. 
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The Green's function (|2.T[) can be obtained under the Lorenz gauge or the Harmonic gauge. However, it can be 
defined under a larger class of gauge conditions which we call by the 'physically reasonable class'. Its definition can 
be understood by considering the residual gauge transformation. By the gauge transformation x a — > x a + £ Q with 
the gauge field of the form 

£ a ( x ) = e- <u < fc - , ""'' +im ^ fe ^ m) ( rj 6»)e 4(few " tr+ ^ 8t8+maj * t * ) , (2.13) 

k,l,rn 

the expansion coefficients h a p( k i m }(r, 9) are transformed, but, the formal expression of the linear metric perturbation 
()2.10p is invariant. 

The main features of the gravitational waveforms can be read off from the formal expression (|2.10[) . For any distant 
observer at r — > oo and a specific angular position, the waveform can be written as 

h(t) = h^j^e-^^-^-^^ 1 ^- 1 ™^-^ . (2.14) 

Because the characteristic frequencies of waves tv r , tog, and are all associated with geodesies, we may conclude 
that the waveform (I2.14| does not reflect dissipative effects occurring on the radiation reaction time scale. This is to 
be expected since the source of the linear metric perturbation moves on a stable geodesic. 



III. SELF-FORCE AND THE ORBITAL EVOLUTION 



Because the gravitational wave part of the linear metric perturbation carries away energy, the orbit of the particle 
deviates away from a geodesic. The effect of the orbital deviation is described by the so-called self-force, that is 
both induced by and acting on the particle itself. Because the metric perturbation induced by the particle diverges 
along the orbit, a regularization prescription is necessary to obtain the finite self- force [5(. Up to leading order in the 
particle's mass, we developed the regularization prescription by the technique of the matched asymptotic expansion. 
The final result is now called the MiSaTaQuWa self- force [f|, and is formally written as 

^-v a := M r(r) = lim ///<* [gW - gW^-](z) , (3.1) 

where the bare term gW represents the full linear metric perturbation induced by the point particle, and the counter 
term gW sm 9- is the singular part of the linear metric perturbation to be subtracted for regularization. / Q [] is a 
derivative operator for the self-force. Because both terms gW and gM sm 9- are divergent along the geodesic, it is 
necessary to evaluate them at a field point x a , not at the geodesic. After the subtraction, gW — gW sm s- becomes 
regular along the geodesic, and we may take the limit x a — > i" to obtain the finite self-force. 

One can use the formal expression (|2.10p to evaluate the bare term gM. The formal expression (|2.10[) is also 
applicable for the counter term gW sm 9- because one can derive the Green's function for the counter term in the same 
form (|2.7|) as we discussed in Ref. 0. After some formal calculations, the self- force can be written as 

M/" - Y.^h e%kXr+UXe ■ (3-2) 

k,l 

It is crucial to observe that the orbital equations (|2.ip and (|2.2p are still true for time dependent £ a , and we will 
use these equations to calculate the orbital evolution due to the self-force. Quantities that evolve due to the self-force 
are denoted with tilde. The evolution equations for £ a are then 

If each of the £ a is expanded in a manner similar to (|3.2p . these equations can be formally integrated to give 

E = E + (e) A + ]T E k ^+^e , L = L + (i) A + ]T L k ^ +il ™ , (3.4) 

C = C + (d) A + J2 C k je lk ^+«x» t (3.5) 

k,l 
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where £g denote the initial values at A = 1 . (£y and £%i are °f or( ier (fi/M) 2 and due to the self-force. 

The orbital evolution can be derived, from (12. ip . (|2.2j) . ()3.4|) and (13. 5|) . Following our notation convention introduced 
above, the inspiralling world line is denoted by x a . We define the r- and 8- motions by 

f = E fke%k * T > ^ = E fa**' ' ( 3 - 6 ) 

k I 



where the expansion coefficients fk and Qi are the same as those in (|2.3I) and (|2.4p . but they are functions of £ a instead 
of £ a . The evolution equations for \ r and xe are given by 

where <i A = (d/dX) acts on £ a of and 0;. The effective frequencies T r and Tg are the same as T r and Tg but are 
functions of £ a instead of£ a . Similarly, we define the t- and 0- motions by 

* = x*+E ir kt ik * r + E > ^ = + E ^ e4fcXr + E $ ei ' xe • ( 3 - 8 ) 

k l k I 

Then, the evolution equations for \t and x<l> are 



d _ 
-Xt 



+J2{ a { d ^ s - *») *? - d ^} e ^ e ' ( 3 - 9 ) 

I 

= f + E{^(rfAXr-T r )^-d A ^}e^ 

+ E { il ( d ^ s - fs ) $ - d ^i} eil * e > ( 3 - 10 ) 



d_ 

dX 



where the effective frequencies T t and are the same as T t and but are functions of £ a instead of £ a . 

We now derive the formal expression of \a by perturbation. Recall that \a — Xa is of order ji/M by definition. 
The formal expressions (13. 7p . (|3.9p and (|3.10p can therefore be written in the form 

^Xa = T„ (0) + (f a )Hp 4 , 1)e «^ + "», (3.11) 

where T a (o) = ^ a \s=s is f° r the background orbital evolution. The formal expression of Xa becomes 

T Q(0 )(A - \ a ) + \ (t Q ) A 2 + T Q(0 , 0) A + Y, X a( k,i)e lkXr+llxa , (3.12) 



Xol ~CK^Uy\' - '* / 1 2 

k,l 

where we set such that \a — Xa when A = 2 . This expression has two key features. One is the linear growth from 

i i^a) ^ 2 ■ This effect comes from the linear growth of £ a since < T Q >= (dT a /d£ a ) < £ a >. This is expected to 

give the dominant phase evolution of the gravitational waveform because of the quadratic growth of the phase. It 
is the reason that the orbit's phase deviates from that of a geodesic on the dephasing time scale (oc 1/ A //I), while 
its frequencies deviate on the radiation reaction time scale (oc l//x). The second feature is a small shift of the time 
averaged frequencies by T Q ( ). The phase evolution due to this effect will accumulate in time 7], however it remains 
small over the radiation reaction time and is therefore not likely to be observable [l|. 



1 We choose ££ , such that J2k,l £ k l e ikXr+ilx e = o when A = 0. 

2 This is achieved by choosing Xa(o,0) such that J^ fc ( Xa(k,l) e Xr "^ e = at A = 0. 
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IV. THE SECOND ORDER METRIC PERTURBATION 



In this section we derive formal solutions of the second order Einstein equation (| 1 . 5|) . By using the same gauge 
condition as we solve (| 1 .4|) . we can formally integrate (|1.5p with the tensor Green's function (|2.7[) . We separate the 
second order metric perturbation into pieces 

[2] [2]N L [2]SF -.s 

9af3 = 9 a p +9 al 3 , (4.1) 

where g a/3 is due to the nonlincarity of the Einstein equation and g aj3 is from the perturbation of the source. 



(II. 5[) then becomes 

G<[l]^[g(2)iVLj = _ G f[2]/ii/r g (l) )g (l)l ( G [l]/x^ g (2)SFj = g^Pl/w . (4 2) 

We first discuss the formal calculation of g^ NL . By the Green's method, we have 

9 l2]NL a p(x) = - [ dx>*G a/3 , v GW[ S W, S W}, (4.3) 

J — oo < r * < oc 

where the domain of radial integration is outside the outer horizon and inside the future infinity; with r* = J dr(r 2 + 
a 2 )/{r 2 ~ 2Mr + a 2 ), this is defined as — oo < r* < oo. Because the Kerr geometry is stationary an axisymmetric, 
the differential operator G^l^h, k] does not include terms that explicitly depend on t and </>. Therefore, the formal 
expression of the source term is 

G<[2]^g(i) ; g(i)] = e -i^ k ,i, m )t+im<f>Ql^ ^ fi g) e '(^.- tr +^' 9 +™/) ; (4.4) 

k,l,m 

where we have used f|2.10D . Using the Green's function of the form (I2.7p . one may obtain the formal expression 

9 [2]NL a p(x) = e-'^t^'.-'^^^yg^^r, ) e i(^^+^* e + m ^^) _ (4. 5 ) 

We see that g^ NL has the same formal expression as the first order term (|2.10[) . This part of the second order metric 
perturbation changes on the dynamical time scale of the orbit which is much shorter than the radiation reaction time 
scale. This means that g^ NL does not include radiation reaction effects, and that g^ NL describes the non-linear 
effect of the Einstein equation. This effect makes a small correction to the wave amplitude of the linear metric 
perturbation, (|2.10|) , and is not likely to be observable. 

Although this conclusion may look reasonable, this argument is not yet well-supported mathematically because 
(|4.3p is actually divergent. In order to evaluate the formal expression of g^ NL , it is necessary to see the possible 
change due to the regularization calculation. The integration of (|4.3[) has two kinds of divergences. One was pointed 
out in Ref. 0], where it was referred to as an ultraviolet divergence. Near the particle, the singular behavior of the 
linear metric perturbation under the harmonic gauge condition is gW oc /i/i?, where R is the spatial distance from 
the particle in the local inertia frame. Because the Einstein tensor is a second order differential operator, we have 
[g(i) ; g(i)j (x fj, 2 /R 4 . The Green's function behaves as oc 1/R. Therefore the spatial integral in (|4.3p diverges as 
l/i? 2 near the orbit. The other kind of divergence, an infrared divergence in the field theory terminology, is due to the 
fact that the graviton is a massless particle. Because the geodesic source of the linear metric perturbation is a stable 
orbit, it can radiate an infinite amount of energy from the infinite past, which is stored in a Cauchy surface. This 
becomes a source for the second order perturbation given by (|4.3[) . At spatial infinity, the linear metric perturbation of 
a specific frequency under the harmonic gauge condition behaves asymptotically as gW ~ /ie~ 4t ^*~ r )/?~ where r is the 
radius in asymptotically flat coordinates. This means that GM^JgW, g' 1 '] ~ /i 2 woj' ' e ~ l ( UJ + UJ )(*- I, )/ r 2 j an ,4 ^ ne S p a tial 
integral of (I4.3[) diverges as ln(r) at infinity. In App. [B] and App. [Cl we introduce the regularization prescription 
for (|4.3|) and show that it does not alter the result. By applying a gauge transformation, the result should apply for 
metric perturbations under broader gauge conditions. 

We next discuss the formal calculation of g[ 2 l SF . In order to identify the stress-energy tensor for the second order 
metric perturbation, we recall the derivation of the self-force by the linear metric perturbation. In Ref. 0, the 
self-force was derived using a matched asymptotic expansion of the black hole metric for the particle and the linear 
metric perturbation away from the orbit. The singular part of the linear metric perturbation around the particle 
(that is, the part of g that diverges as n/R) is correctly matched to the corresponding singular behavior of the black 
hole metric. In order to match the regular part of the linear metric perturbation with the corresponding part of the 
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black hole metric, it is necessary to account for the motion of the black hole with respect to the local inertia frame of 
the background metric. The result is the regular equation of motion. This suggests that the stress-energy tensor to 
the second order perturbation can be written as 

J V^m 

where x a and v a are the position and 4-velocity of the orbit derived in Sec. IIIII 
The metric perturbation is then 



J + ffLT F = 8 ^ / d\ (^j G a p , M „ (x, 5(A)) ~v»v v . (4.7) 

Recall our definitions (|3 . 6|) and (|3.8[) for the orbital coordinates. Here, we use the same formal expressions (|2.3[) - (|2.6[) 
for the geodesic, but we replace £ a and Xa by £ a and Xa- Thus we can rewrite (|4.7p by replacing £ a and Xa of (|A3j) 
with £ a and Xa 

9ap+9 l a/ F = 87r/i J dX e-^ t+4m ^^' i ' m (r,6»;£)e ,; ^- im ^- l ' c ^-^ ; ( 48 ) 

Lj,k,l,m 

where h^' l ' m is defined in the same way as h^'p' ' m > but with £ a replaced by £ a . 

We now discuss the effect of the self-force on (I4.8[) by perturbation. There are two kinds of effects. One is from the 
dependence of the expansion coefficients fi^' l ' m f £ a , and the other is from the phase function Xa in the exponentials. 

We first look at the effect of the expansion coefficients hf^'p ' ' m . Using (|3.4p and (|3.5p , the expansion coefficients 
become 

h»J> l > m (r,e;£) -> Kf- m (r,8;£ Q ) + A_h^ l - m (r,6;£ Q ) A, (4.9) 

where we have ignored the oscillating part of £ a . (The metric perturbation due to the oscillating part of £ a can 
be written in the form (|2.10p , and does not show any new characteristic feature.) One can see that the expansion 
coefficients grow linearly with A. This is because of the linear growth of £ a caused by the self-force. This leads to the 
linearly growing/decaying feature of the wave amplitude through the £ "-dependence of h a p' l ' m) ' ' ' as 3 

^ n a^., \ " „-iwtk i m \t+imd> ( i^(k,i, m )M-l,ra , s ;w (M _ m) ,fc,Z 



T 



t(0) 



87T/X J2 e-^^ t+m * (h^ l - m) ' k,l ' m (r, 0) + fc$' , ' m "*' ,,m (r, 0)tj e ^e+l^m^) ^ (4 n) 



k,l,m 



where h^ 1 ^' " = {dh^f l ' m / d£ a ) (£ a ^ /T t . This formal expression for the waveform has a new feature, that is, 
the wave amplitude changes on the radiation reaction timescale. 

We next see the effect of Xa in the exponential of (I4.8[) . Using (|3 . 1 2[) this exponential has a quadratically growing 
feature 



JuiXt -irriX4,-ikXr -USce giwf t(0) A-irnT^o) A-ifcT r(Q ) \-UY e(0) A g-«"T t(o) A'+imT +( o) A^+ifcT^o) A''+i2T e(0) A 8 

X (l + (iu (t t ) - im (f \ - ik (t r ) - il (to)) , (4-12) 



where we have ignored the oscillating part of Xa (Again, the metric perturbation due to the oscillating part of Xa 
can be written in the form (|2.10[) . and does not show any new characteristic feature.) and we use the renormalized 



The actual procedure of integration is 

J dwe-" 1 ^ J dt't'e^-^' = J dwe-^i-i)^- J dt'e 1 ^-^' = (2tt) J due'^ {-i)-^-5{u - a) ) 

= (2tt) f dwStui - ui )i — e~ iwt = (27r)te- <Wo * . (4.10) 
J did 
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frequencies T Q ( ) = T Q ( ) + T Q ( ,o)- We obtain the formal expression as 



4 



where we used 



W( fc ,i, m ) = rn^- + k-^ + , (4.15) 

J-t(o) J-tCo) J-tco) 





T0( O ) (T t ^) 


1 t(0) 


1 t(0) 


■<*•) 




-f2 

1 t(o) 


1 t(0) 




W(fe,z irra ) = rn 



\ 1 *(o) i t(o) y 

From the time-dependence of the exponent, one can see two key features of the wave phase. First, the quadratic 
growth of the wave phase on the dephasing time due to the linear growth of the fundamental frequencies. Second, 
the linear growth of the wave phase over the radiation reaction time due to the renormalization of the fundamental 
frequencies. Both features have been seen from the orbital evolution due to the self-force as we argued in Sec. Mil 
The results of (|4.ll|) and (|4. 14[) suggest that the gravitational field might be well-approximated by the form 



k,l 



h %' m ( r ' d > ^)) e- lk ^ {t) - ll ^ {t) ' lm ^ l) , (4.17) 



where £ and = r, 8, cf>) are functions of the Boycr-Lindquist coordinate time t and the radiation reaction effect is 
included in the waveform through those functions. 



V. CONCLUSION 



In this paper, we discuss the modulation of the gravitational waveform due to gravitational radiation reaction. Since 
the linear metric perturbation is induced only by the background geodesic, it is necessary to calculate the second order 
metric perturbation in order to see the radiation reaction effect. A more rigorous calculation of the second order metric 
perturbation is not yet available. For this reason, we have considered only a qualitative study, extending the technique 
of formal calculation used to obtain the linear metric perturbation in Ref. Q ■ The advantage of this technique is 
that one can grasp some key features of the gravitational waveform without a complicated calculation. 

There are two kinds of source terms for the second order metric perturbation. One is the non-linear term of the 
linear metric perturbation G^^lg^, g^]. Because the formal calculation of this term is triperiodic, the part of 
the second order metric perturbation (|4.5[) is also triperiodic. The waveform induced by this part has exactly the 
same spectral form as that of the linear metric perturbation. This part simply changes the amplitude of the waves, 
therefore, describing the correction to the wave propagation due to the non-linearity of the Einstein equation. The 
other kind of source term is the point source due to the orbital deviation from the background geodesic, T^ 2 ^ . The 
part of the the second order metric perturbation induced by this source has two new features; the linear growth of 
the wave amplitude over the radiation reaction time scale [See (I4.HI) .] and the quadratic growth of the wave phase 
over the dephasing time scale [See (|4.l4l) .]. From the derivations of (|4. 1 1|) and (|4.l4p . it should be clear that these 
features reflect the secular effect of the orbital evolution due to the self-force discussed in Sec. Mil 



The actual procedure of integration is 

J due- iut J dt't' 2 ^^"^' = J dwe-^i-l)-^ J dt'S"-" ^' = (2tt) J due- iut (-1) 5(u> - u ) 

= (2tt) / dw5{u> -£j )(-l) — -e'^ 1 = {2-K)t 2 e~ l ^ t . (4.13) 
J dw 2 
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These results suggest that the effect of radiation reaction on the waveform are included in templates of the form 
(|4.17p , which is the linear metric perturbation with the geodesic constants replaced by those which evolve adiabatically 
due to the self-force. This result agrees with the procedure of the adiabatic approximation. However, strictly speaking 
the perturbation scheme is valid over a dephasing time and over that time scale, only the feature of the quadratic 
growth of the wave phase can be identified as a secular effect. Beyond this time scale, it is necessary to account for 
the secular effect of the self-force due to the second order metric perturbation. One can only discuss those linearly 
growing features in a gauge invariant way after calculating the third order metric perturbation. 

Because the perturbation scheme we use is limited to the dephasing time scale, it is not an appropriate method 
for calculating gravitational waveforms. In order to solve this problem, we propose an adiabatic expansion that is 
a systematic perturbation of a held theory coupled to a particle Because this expansion recovers the adiabatic 
approximation at leading order, our result here supports the validity of this new expansion method. The application 
of this technique to the adiabatic expansion shall be discussed elsewhere. 
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APPENDIX A: FORMAL EXPRESSION OF THE LINEAR METRIC PERTURBATION 

In this appendix, we review the formal calculation of the linear metric perturbation given as (|2.9|) . With the ansatz 
(|2.7p . we decompose the Green's function as 

G a f3^{x, x') = ^9^( r ' r/ ^' e ')e- iw(i - t ' )+im( ^' ) . (Al) 

uj.in 

We plug this into (|2.9p and we have 

g%{x) = STr/i f dX f^) ^^(r 1 f;M>' t /e-( i - i )+ im »-« . (A2) 

At this point, we recall the periodic feature of the orbit (|2.3l) - (|2.6p . It is easy to see that dr/dX and 9a'aliu( r '^^^) 
can be expanded as a discrete Fourier series in e~ lkXr ~ dxe . From (|2.ip and ()2.2j) . v^v v can be similarly expanded. 
As for e tuJt ~ m 4> ( one ma y expand by the same discrete Fourier series except the factor e l " x * _lmx *. In summary, one 
can formally rewrite (|A2[) as 

By integrating over A, we obtain 

g m( x ) = Y 8n ^ J2 e-^^ t+im ^h^ k - l - m) ' k ' l ' m (A4) 

which is equal to (|2.10[) . 

APPENDIX B: REGULARIZATION OF THE ULTRAVIOLET DIVERGENCE 

We use the regularization prescription for the ultraviolet divergence proposed in Ref. By the consistency of 
the matched asymptotic expansion, we know that the divergence behavior of the second order metric perturbation 
becomes g}- 2 ^ NL cx /i 2 /R 2 , with respect to the local inertial coordinates in the neighborhood of the geodesic. (R is the 
spatial distance from the particle in the local inertia frame.) The idea in Ref. [8| is to subtract this divergence by 
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using the quadratic combination of the scalar field induced by a point particle and calculate the remaining regular 
part by the Green's method. We use the scalar field that satisfies 

under the retarded boundary condition, where x^(t) is the same geodesic used in (|2.8j) for the linear metric pertur- 
bation. The resulting scalar field has the divergence behavior $ oc 1/R near the geodesic. With appropriate tensors 
^ajs an d which are regular around the orbit, one may construct 

such that g[ 2 ] JYi and g( s ^ 2 ^ NL have the same singular behavior near the geodesic. We define the remaining part by 



which satisfies 



(R)[2]NL _ [2]NL (S)[2]NL , m \ 
— 9 a fj 9 a /3 ' \ a; >) 



Gf[i]/^r g (-R)(2)wii = _ G [iWr g (S)(2)iv-ii _ G [2]^r g (i) )g (i)i _ m 4 \ 
Because g^g^ NL is regular, the RHS of (|B4[) must be regular near the geodesic, and g( R )( 2 ) NL can be written as 

g mWL {x) = f d*x>G a ^ v f G mns {SmNL ]+G^S {1 \s^) ■ (B5) 

J— oc<r*<oo 



In order to see the formal expression for (|B5|) , it is necessary to know the formal expression for (|B2[) . The expression 
for $ can be obtained in a similar manner as (|2.10|) . Because the scalar Green's function of (|B1[) can be written in 
the form 

G(x, x') = G(t -t',<f>- (/)'; r, r', 9, 9') , (B6) 
similarly to (|2.7p , the formal expression of $ becomes 

= ^ e- l "< fc - i -™)* +lm0 $ (fc . i!m) (r,6';£ a )e l(fcw '' tr+i ' i ' 8 * £ ' +m ^ t * ) . (B7) 

k,l,m 

The structure of the singularity depends only on the local background geometry, and the components of and Jfcj^g 
can be constructed only from the background curvature along the orbit. Therefore, we have 

k,l,rn 

for n = 1, 2. In summary, the RHS of (|B4[) has the same formal structure as (|4.4|) . 



APPENDIX C: REGULARIZATION OF THE INFRARED DIVERGENCE 

Though we can remove the ultraviolet divergence by calculating (|B5[) instead of calculating (|4.3|) . it does not 
guarantee that we can also remove the infrared divergence, i.e. the integration over the spatial volume of (|B5|) 
becomes divergent at the large spatial radius and at the black hole horizon. 

A formal method of avoiding this divergence is to introduce the radial regulator T to (|B5|) as 



gif 2]NL (x;T) = [ dx'^G^ ( G ^[^ NL ] + G^[^\ g« 

J-r<r*<r v 



(CI) 



By definition, g (. R )l 2 \ NL (x; T) satisfies (TB4|) for -T < r* < T, and it diverges as g ( R )\ 2 \ NL (x; T) cx ln(T), when T -> oo. 

The counter term to cancel this divergence g ^ c ^ NL (x; T) must satisfy the vacuum linearized Einstein equation 
so that g(- R )[ 2 ] JVi + g( c )l 2 ] NL still satisfies (|B4[) . Because g^ NL is intrinsically divergent along the orbit, one may 
consider a point source moving along the geodesic which induces g^ c ^ NL (x; T) as 

T^ NL ^(x;T) = ln(I> 2 J drP^ 5 -^^-, (C2) 
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where x a (r) is the same geodesic as in (12.8[) . We have the counter term 

g (C)[2]NL {x . T) = ^ J dx >4 Ga ^ {X]X i )T mNL^ {x . T) 

= 87r/i 2 ln(r) ( dTi^G af} ^{x,x). (C3) 



Here fl- 2 ^" is the tensor derivative operator defined along the orbit such that (g(- R )[ 2 1 JVL + g( c )[ 2 ] NL ^ (x;T) becomes 
finite everywhere except along the geodesic when r — * oo. 

We discuss the formal expression of limr^oo (g( R ^ NL + g( c ')[ 2 ] Jvi ^ [x; T) following this regularization prescription. 
For a finite T, as we may see from ()Cip . g( R '^ NL (x;T) has the same formal expression as (|4.5|) because the source 
term of (|4.5j) has the same formal expression as (|4.4|) . By g^°^ 2 ' NL (x;T), we simplify subtract the linear divergent 
piece of g( R ^ NL (x;T), i.e. we have 

„(R)[2]N L ( rr . v i\ 

g (CMNL {x . T)= _ ln(r) ^ 8 (C4) 

and this operation does not alter the formal expression. Thus, by the infrared regularization, the part of the second 
order metric perturbation due to the non-linear term of (|1.5[) has the formal expression of (|4.5p . 



[1] Y. Mino, Prog. Theor. Phys. 115 (2006) 43 ( [gFqc /0601019), Y. Mino, Prog.Theor.Phys. 113 (2005) 733 | |ir^qc /0506003 ), 

Y Mino, R. Price, in preparation. 
[2] Y. Mino, Phys. Rev. D67 (2003) 084027 ( gr-qc/0302075] ). 
[3] S. Drasco, in preparation. 

[4] Recent development of the perturbative approach to the gravitational radation reaction is reviewed in Class. Quant. Grav. 
22 (2005) S543-S859. 

[5] Y. Mino, M. Sasaki, T. Tanaka, Phys. Re v. D55 (1997) 3457, E. Poisson, Living Rev. Relativity 7 (2004), 6 ([Online article] 
|http: / /www.livingreviews.org/lrr-2004-6 ) , E. Poisson, Class . Quant. Grav. 21, R153 (2004). 

[6] L. Barack, A. Ori, Phys. Rev. D 61, 061502(R) (2000) ( gr-qc/99120To|, L . Barack, Phys. Rev. D 62, 084027 (2000) 
( |gr-qc/0005042| , L. Barack, Phys. Rev. D 64, 084021 (2001) (gr-qc/0105040 ), L. Barack, Y. Mino, H. Nakano, A. Ori, 
M. Sasaki, Phys. Rev. Lett. 88 (2002) 091101, L. Barack, A. Ori, Phys. R ev. D 66, 084022 (2002) (gr-qc/020409), Y. Mino, 
H. Nakano, M. Sasaki, Prog. Theore. Phys. (2003)108 ( |gr-qc/0111074| , S. Detweiler B. F. W hiting, Phys. Rev. D 67, 
024025 (2003) ( |gr-qc/020 2086), L. Barack, A. Ori Phys. Rev. D 67, 024029 (2003) ( |gr-qc/02090"72| ), L. Barack, A. Ori, 
Phys. R ev. Lett. 90, 111 101 (2003) dgr-qc/0212lMt, S . Detweiler, E. Messaritaki, B. F. Whiting, Phys. Rev. D 67, 104016 
(2003), (gr-qc/0205079), D. H. Kim, gr-qc/0402014 ( W. Hikida, H. Nakano, M. Sasaki, Class. Quant. Grav. 22, S75 3 
(2005) ( |gr-qc/04m"50X W. G. Anderson, E. E. Flanagan, A. C. Otte will, Phys.Rev D71 (2005) 024036 | |gr-qc/04120"09| , 
W. G. Anderso n, A. G. Wisem an, Class. Quant. Grav. 22 (2005) S783 ( |gr-qc/0506136| , L. Barack, N. Sago, Phys. Rev. D75 
(2007) 064021 ( |gr-qc/0701069j ), R. Haas, Phys. Rev. D 75, (2007) 124011 (gr-qc/0704.0797). 

[7] A. Pound, E. Poisson, gr -qc70708.3037, A. Pound, E. Poisson, gr-qc/0708.3033, A. Pound, E. Poisson, B. G. Nickel, Phys. 
Rev. D72 (2005) 124001 ( |gr-qc/0509i22| |. 

[8] E. Rosenthal, Phys.Rev. D74 (2006) 084018( |gFqc /0609069 ) , E. R osenthal, Phys.Rev. D73 (2006) 044034 Hgr-qc/060 2066 ), 
E. Rosenthal, Phys.Rev. D72 (2005) 121503 (gr-qc/0508050 ), E. Rosenthal, Class.Quant. Grav. 22 (2005) S859 
( |gr-qc/0501046p . 



